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Abstract. We present observations of nine radio pulsars using the Heinrich-Hertz- Telescope at A 0.87 mm and the 
IRAM 30-m telescope at A 1.2 mm in search for a cold dust around these sources. Five of the program pulsars 
have been observed for the first time at the mm-wavelengths. The results are consistent with the absence of 
circumpulsar disks that would be massive enough (> O.OIMq) to support planet formation according to the 
scenarios envisioned for solar-type stars, but they do not exclude lower mass (< 10 — lOOM^) disks for a wide 
range of grain sizes. These conclusions confirm the previously published results and, together with the current lack 
of further detections of pulsar planets, they suggest that planet formation around neutron stars is not a common 
phenomenon. 
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1. Introduction 

The existence of planets around one of the millisecond pul- 
sars, PSR B1257-I-12 (Wolszczan & Frail 1992; Wolszczan 
1994), has raised a possibility that at least some pulsars 
may, like normal stars, be accompanied by protoplane- 
tary or debris disks. Most theories of planet formation 
around neutron stars assume a creation of some sort of a 
protoplanetary disk out of the material supplied by the 
pulsar's binary companion or, possibly, by the fallback 
of supernova ejecta (e.g. Podsiadlowski 1993; Phinney & 
Hansen 1993). Of course, properties and composition of 
such disks may be significantly different from those of the 
disks commonly observed around young normal stars (e.g. 
Boss 2003). In addition, the initial conditions for planet 
formation in an expanding disk must depend on its abil- 
ity to protect itself from high-energy photons generated by 
early accretion onto the pulsar and from a wind of ultra- 
relativistic particles that carries most of the pulsar's spin- 
down energy (Phinney & Hansen 1993; Miller & Hamilton 
2001). Although the spindown luminosities of pulsars are, 
in principle, more than sufhcient to heat the dust grains, 
the efficiency of this process depends on a degree of cou- 
pling of the pulsar wind to dust, which is difficult to pre- 
dict in absence of a sufficient observational evidence. 

Since the discovery of the PSR B1257+12 planets, 
there have been several attempts to detect the hypo- 
thetical circumpulsar disks with both the space- and 
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the ground-based telescopes. Upper flux limits have been 
derived for a number of sources at wavelengths rang- 
ing from 10 /im to 3 mm (van Buren & Terebey 1993 
(IRAF); Zuckerman 1993 (IRTF); Phillips & Chandler 
1994 (JCMT, OVRO); Foster & Fischer 1996 (IRTF); 
Greaves & Holland 2000 (JCMT); Koch-Miramond et al. 
2002 (ISO, IRAF/Scanpi); Lazio & Fischer 2004 (ISO)). 
As shown by these authors, the existing infrared and mm- 
wave flux limits for PSR B1257+12 and other pulsars do 
rule out massive, > O.OIM© disks similar to those thought 
to give rise to planets around normal stars (Boss 2003). 
However, these limits do not contradict a possibility that 
some pulsars may be accompanied by much less massive 
disks ranging from a fraction of the asteroid belt mass 
to a few hundred M©. In fact, the three PSR B1257+12 
planets of a total mass of ~8 M,^ appear to have formed 
from a protoplanetary disk, as suggested by the observed 
coplanarity of their orbits (Konacki & Wolszczan 2003). 
Along these lines, Miller & Hamilton (2001) have argued 
that the lower limit to a disk mass necessary to shield 
the forming PSR B1257+12 planets against the particle 
and the photon emission of the pulsar would be around 
^10^^ g. Furthermore, Hansen (2002) has considered PSR 
B1257-fl2 planet formation scenarios in which the ex- 
panded, planet building part of the disk would require 
a mass comparable to that of the existing planets. 

In this paper, we present observations of nine pul- 
sars with the Heinrich-Hertz- Telescope at A 0.87 mm and 
the IRAM 30-m telescope at A 1.2 mm. Our measurements 
were designed to either detect or to set upper limits to an 
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Table 1. Continuum flux limits for nine pulsars at A0.87mm and A 1.2mm 



PSR 


P 


P 


d 


r 


So. 87 


Sl.2 


Remarks 




(ms) 


(10-^«ss-i) 


(pc) 


(Myr) 


(mJy) 


(mJy) 




J0108-1431 


807.6 


77.0 


130 


166 


14.5 


— 


solitary PSR 


B0950+08 


253.1 


229.8 


160 


17.5 


14.6 





solitary PSR 
MSP- WD binary 


J1012+5307 


5.3 


0.017 


520 


4860 




2.7 


J1024-0719 


5.2 


0.019 


350 


4410 


— 


1.8 


solitary MSP 


B1257+12 


6.2 


0.114 


620 


863 


7.6 


2.0 


MSP, 3 planets 


B1855+09 


5.4 


0.018 


700 


4760 


12.5 


2.0 


MSP- WD binary 


J2124-3358 


4.9 


0.021 


250 


3800 


56.0 


— 


solitary MSP 


J2307+2225 


535.8 


8.7 


380 


976 


12.0 


— 


solitary PSR 


J2322+2057 


4.8 


0.010 


780 


7850 


10.2 


— 


solitary MSP 



Notes.- (1) Pulsar distances, d, were derived from the Taylor & Cordes (1993) model of the Galactic electron density distribution. 
(2) The characteristic age of a pulsar is defined as r = P/2P, where P is the pulsar period and P is the period derivative. (3) 
The flux values, So. 87 and Si. 2, are 2a limits. 



emission from possible cool protoplanetary or debris disks 
around these objects. Five of the program pulsars have 
been observed at these wavelengths for the first time. We 
describe the instrumental setup and present flux density 
measurements of the nine sources in Sect. 2. In Sect. 3, 
these and other results are further discussed in the context 
of the recent observational and theoretical developments. 

2. Observations and data analysis 

The sample of pulsars selected for our survey consisted 
of nine nearby sources within the distance limit of < 1 
kpc. Alongside with the planets pulsar, PSR B1257+12, 
we have included three other solitary millisecond pulsars 
(MSPs) on the assumption that they would be the most 
probable candidates for a detectable dust emission. As bi- 
nary MSPs must have undergone long periods of accre- 
tion from their companions and could be left with a resid- 
ual circumstellar material suitable for planet formation 
(Phinney & Hansen 1993), two MSP- WD (white dwarf) 
binaries have been also included. Finally, for complete- 
ness, we have supplemented our list with the three nearby 
normal pulsars that would be observable at the telescope 
sites. 

The continuum observations at A 0.87 mm were car- 
ried out at the Heinrich-Hcrtz- Telescope (Baars et al. 
1999) on Mt. Graham, Arizona, during two observing 
sessions in January 2001 and December 2002. We used 
a sensitive 19-channel bolometer array developed by E. 
Kreysa and collaborators at the Max-Planck-Institut fiir 
Radioastronomic, Bonn. The centre frequency of the 
bolometer was about 345 GHz, with a bandwidth 70 GHz. 
The beam size on the sky was ^^23". Using the central 
channel of the bolometer, we performed standard on-off 
measurements with the telescope secondary chopping in 
azimuth by 40". In order to calculate the atmospheric 
zenith opacity, we made skydip observations every 45 min. 
The resulting opacities at the observing frequency varied 
between 0.3 and 0.9 during the two observing sessions. 



For calibration purposes we performed mapping and on- 
off measurements of various planets (mainly Mars, Saturn 
and Uranus), resulting in conversion factors of 7.7 — 9.5 
mJy per count. 

The observations at A 1.2 mm were carried out in 
February and March 2003 at the IRAM 30-m telescope on 
Pico Veleta (Spain). We used the 117-channel Max-Planck 
Millimeter Bolometer (MAMBO-2) array (Kreysa et al. 
1999), which had a half-power spectral bandpass from 210 
to 290 GHz with an effective frequency of 250 GHz. The 
beam size on the sky was 10.7 arcsec. The sources were 
observed with a single channel using the standard on-off 
mode with the telescope secondary chopping in azimuth 
by 40" at a rate of 2 Hz. The sky opacity at the observing 
frequency was monitored with skydips and varied between 
0.1 and 0.3 during the seven days of the observing run. 
For flux calibration we observed a number of calibration 
sources, resulting in conversion factors of 32 — 37 counts 
per mJy and an estimated absolute flux uncertainty of 
15%. 

The data obtained at the two telescopes were con- 
verted into a standard format and analysed using the NIC 
software package. Correlated noise was subtracted from 
the reference channel using the weighted average signals 
from the surrounding channels. Anomalous signals were 
clipped above the 5cr level. 

We did not detect dust emission towards any of the 
pulsars observed. Table 1 lists the 2a flux limits obtained 
at the two observing frequencies. At A 0.87mm we reached 
limits of 7 — 14 mJy for most of the sources, which is 
slightly more than the limits of the survey by Greaves & 
Holland (2000). At A 1.2 mm 2cr flux limits of 1.8-2.7 mJy 
were obtained. 

3. Discussion 

A lack of detections and the derived flux limits at 
A 0.87 mm and A 1.2 mm for the nine pulsars discussed in 
this paper (Table 1) place useful constraints on the exis- 
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Fig. 1. Observational limits and theoretical models for the emission from dust in the PSR B1257+12 planetary 
system. Symbols mark the flux limits set by observations discussed in the text (o) and by published measurements 
with the following instruments and surveys: IRTF(*,o), IRAS/Scanpi (•), JCMT/OVRO (A, x), and 2MASS (*). 
The solid curves trace the loci of a maximum flux for theoretical disk emission models calculated for four disk masses 
as a function of the grain size ranging from 0.01 to 1000 /im. Dashed lines correspond to the Spitzcr Space Telescope 
3cr sensitivity limits estimated for the IRAC and the MIPS detectors. 



tence of cold, dusty disks around these sources. To illus- 
trate the implications of this result for a circumpulsar disk 
observability in the framework of a simple model, we use 
the approach of Foster & Fischer (1996) by assuming that 
a fraction / of the pulsar spindown luminosity E heats a 
disk of mass nid consisting of dust grains of size a to a 
temperature T, and calculate the expected infrared flux 
as a function of disk parameters and pulsar distance fol- 
lowing the formalism developed by these authors. Because 
PSR B1257+12 has planets that have originated in a disk 
(Konacki & Wolszczan 2003), it is the best sampled pul- 
sar in the sub-mm/mm range, and the published disk mass 
limits are not dramatically different from one pulsar to an- 
other, we can use this object as a representative example 
to make a general comparison of model calculations with 
the data. 



In Fig. n four sets of models for disk masses m^ = 
10""*, 10"^ 10", and WMq, are compared with the flux 
limits for PSR B1257-I-12 obtained in this work and the 
ones published elsewhere. To ensure compatibility with 
other results, we have assumed fc=0.01 (Foster & Fischer 
1996), 0.01 /xm < a < 1000 ^m (Koch-Miramond et al. 
2002), the standard spindown luminosity of the pulsar, 
E = 2 X 10^* erg s~^, and the pulsar distance of 620 pc 
from the Taylor & Cordes (1993) model of the Galactic 
electron density distribution. In the modelling process, for 
the assumed disk mass and grain size ranges, the grain 
temperature varied over the range of 8 K < T < 500 
K, from the largest to the smallest values of rud and a, 
respectively. This is practically within the range between 
the typical temperature of the cold interstellar dust (10 K) 
and the sublimation temperature of silicate grains (1500 
K) (e.g. Koch-Miramond et al. 2002). 
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A comparison of the existing observations with a sim- 
ple dust model shown in Fig. ^ clearly demonstrates the 
absence of a disk around PSR B1257+12 with the mass 
> O.OIMq that appears to be necessary to support planet 
formation around solar- type stars (Boss 2003). This con- 
clusion is also valid for our data on other pulsars, and 
it is generally consistent with the limits on dusty disk 
masses obtained from previous observations of a number 
of pulsars in the AlO/xm - A 3 mm range (e.g. Greaves 
& Holland 2000; Koch-Miramond et al. 2002). Of course, 
one could alter this outcome by varying the heating ef- 
ficiency of the disk by the pulsar wind, or changing the 
disk model assumptions. For example, an upper limit in 
excess of O.OIMq can be obtained for the mass of a disk 
around PSR B1257-f 12 by extrapolating from models of 
the T Tauri disks (Phillips & Chandler 1994). Similarly, 
a dusty disk model used by Jura (2003) to explain an in- 
frared excess in the white dwarf G29-38 can, in principle, 
be modified to describe the illumination of a disk by non- 
thermal radiation from a neutron star (M. Jura, personal 
communication). However, these results cannot be con- 
clusively compared with the ones based on more generic 
models (e.g. Foster & Fischer 1996; Koch-Miramond et al. 
2002) without a much better understanding of the physics 
of hypothetical circumpulsar disks. 

As stated earlier and evident from Fig. ^ the data do 
not exclude a possibility that dusty disks ranging from 
a fraction of the mass of the solar system asteroid belt 
up to possibly ^^lOOM^ may exist around pulsars over 
a wide range of temperatures and grain sizes. Conditions 
for the existence of such disks have been recently exam- 
ined by Miller & Hamilton (2001), who specifically predict 
that a leftover material cannot exist around pulsars with- 
out planets. This is because such pulsars have evidently 
not managed to develop protoplanetary disks with a suf- 
ficient mass {rrid > lO^^g for PSR B1257-hl2) to protect 
themselves against evaporation by the pulsar wind and 
accretion flux that would prevent planet formation. On 
the other hand, pulsars with planets, like PSR B1257-I-12, 
could still have observable debris disks left over from the 
final stages of a protoplanetary disk evolution (e.g. Hansen 
2002). In principle, dramatically more sensitive observa- 
tions of PSR B1257-f-12 and other pulsars with the Spitzer 
Space Telescope will have a potential to either reveal such 
a circumpulsar debris material or to place very meaningful 
constraints on its existence (Fig.^. 

Our data support a general conclusion drawn from the 
previous negative results of searches for dust emission from 
the vicinity of pulsars that planet formation around neu- 
tron stars occurs infrequently and that it follows an evolu- 
tionary path that must generally be different from the one 
envisioned for planets around normal stars. This conclu- 
sion particularly refers to possible sources of a circumpul- 
sar matter, as well as to the problem of a disk formation 
and retention around a neutron star (Podsiadlowski 1993; 
Phinney & Hansen 1993; Miller & Hamilton 2001). It is 
also consistent with a crude statistic of < 10% derived 
from the fact that PSR B1257-I-12, one out of the nine 



known solitary millisecond pulsars, has planets around it 
(Konacki & Wolszczan 2003). Similarly, it is not surpris- 
ing that no massive dusty disks have been found around 
relatively few normal, young pulsars observed so far, as 
the formation and survival of planets around such objects 
may be particularly difficult (Thorsett & Dewey 1993). 

In the absence of new detections of planets around 
millisecond pulsars (the Jupiter-mass planet in the M4 
globular cluster, recently confirmed by Sigurdsson et al. 
(2003), was probably created around a normal star) fur- 
ther searches for dust around these objects, covering all 
the physically plausible parameter space at a possibly high 
sensitivity level, is a logical way to meaningfully constrain 
the models of a creation and evolution of pulsar proto- 
planetary and debris disks. Given the existing correlation 
between stellar metallicity and the frequency of planets 
around normal stars (Santos et al. 2003; Fischer et al. 
2004), it is likely that the protoplanetary disks around pul- 
sars, possibly made of a highly evolved stellar material, do 
not have to be very massive to support an efficient planet 
formation (see also Lazio & Fischer 2004). Furthermore, 
one would also have to include a distinct possibility that 
the "dusty plasma" environment of a circumpulsar disk 
(e.g. Mendis & Rosenberg 1994), driven by the relativistic 
pulsar wind and the magnetic fields carried with it, can- 
not be satisfactorily approximated by the standard pro- 
toplanetary or debris disk models based on theories and 
observations of the formation regions of solar-type stars 
(e.g. Boss 2003). In the near future, the Spitzer Space 
Telescope, with its factor of ~10^ — 10^ improvement in 
sensitivity compared to the instruments previously used 
to search for dust emission from around pulsars (Fig. P) , 
is an obvious choice for further exploration of the physics 
of neutron star planetary systems. 
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